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Abstract 

Azimuthal angle dependence of the pion source radii was measured ap¬ 
plying the event shape selection at the PHENIX experiment. The measured 
final source eccentricity is found to be enhanced when selecting events with 
higher magnitude of the second-order flow vector, as well as the elliptic flow 
coefficient V2. The spatial twist of the particle-emitting source was also ex¬ 
plored using a transport model. Results indicate a possible twisted source in 
the final state due to the initial longitudinal fluctuations. 


1 Introduction 

Higher-order flow coefficients v n are useful observables to constrain the properties 
of the quark-gluon plasma, such as a shear viscosity over entropy density ratio, in 
heavy ion collisions mm- The v n are found to largely fluctuate in each event even 
if the collision centrality is fixed [3j, which is due to the initial spatial fluctuation 
of participating nucleons. To control such event-by-event initial fluctuations, the 
event shape engineering was suggested (4]- The event shape engineering can be 
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performed by selecting the magnitude of the event flow vector Q 2 . 


Q'2 = \Jvlx + <l2,y/v r M, (!) 

92 ,x = ycos(20j), q 2 . y = y^WjSin(2^), 

where M is the particle multiplicity, fa is the azimuthal angle of particles, and Wi is 
a weight for particle i. The magnitude of Q 2 is proportional to the strength of event- 
by-event v 2 and also reflects the resolution of the second order event plane. This 
technique could allow us to test the effect of the initial geometry on the final state 
distribution of emitted particles, which help to understand the medium response 
through the system evolution. 

The event shape engineering focuses on the fluctuations in the transverse plane. 
It is naturally considered that there would be the fluctuations not only in the trans¬ 
verse plane but also in the longitudinal direction. The presence of such fluctuations 
could cause a twisted source along the longitudinal direction |5j- The number of 
participants going to the forward and backward directions would be different, and 
therefore participant planes might also be different. As a result, the event plane 
angle in the final state could be different between the forward and backward an¬ 
gles. Thus the initial spatial twist may survive in the final state as well as a twisted 
flow [8] [9]. 

In this proceedings, we present results on HBT measurements using charged 
pions and applying the event shape engineering technique for Au+Au collisions at 
fas NN =200 GeV recorded with the PHENIX experiment. Also, we examine the 
possibility of a spatially twisted source in the final state using HBT interferometry 
in AMPT model. 


2 The event shape engineering at the PHENIX 

The flow vector was determined by the Reaction Plane Detector (RXN, 1 < M < 
2.8). In our analysis, the w, in Eq. ([ 2 ]) reflects the multiplicity or energy measured 
in a segment i of the RXN and M is 'Yh VJ i- The tracking for charged particles was 
performed by the Drift Chamber and the Pad Chamber, and the particle identifica¬ 
tion was done by the electromagnetic calorimeter, where they have the acceptance 
of |?y| < 0.35 and \fa < n/2 in the west and east central arms. 

The measured Q 2 distributions were fitted with the Bessel-Gaussian function [4], 
and then the events were classified with the magnitude of Q 2 . Figure |T] shows 
charged hadron v 2 for the higher 20% and the lower 30% Q 2 events. Results 
without Q 2 selection are also plotted for the comparison. The effect of Q 2 selection 
on the v 2 is clearly observed, that is, the higher (lower) Q 2 enhances (decreases) 
the strength of v 2 compared to the v 2 without the Q 2 selection. 

Then the Q 2 selection was applied to the HBT measurements using charged 
pion pairs. In the HBT analysis, pion pairs were analyzed with the out-side-long 
parameterization d QA] in the longitudinally co-moving system. The effect of 
the event plane resolution was also corrected for both cases with and without Q 2 
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Figure 1: Charged hadron t; 2 as a function of pt for four centrally bins with and 
without the Q 2 selection in Au+Au 200 GeV collisions. The lower and higher 
Q 2 events were selected. 


selection. Figure |2^left) shows the extracted pion HBT radii, R% and R°, as a 
function of azimuthal pair angle (f> relative to the second-order event plane 4' 2 . 
Results show that the higher Q 2 selection increases the oscillation strength compared 
to the case without Q 2 selection. 

These oscillations of HBT radii are supposed to be sensitive to the final source 
eccentricity at freeze-out, eg na i. Blast-wave studies suggest that the quantity of 
2R^ 2 /R s ,o would be a good probe of £g na i in the limit of fcx = 0, where kj- denotes 
a mean transverse momentum of pairs. The oscillation amplitudes of R% and in a 
form of the final eccentricity are plotted as a function of the number of participants 
calculated by Glauber model in Fig. fright). The higher Q 2 selection enhances the 
measured £g na i as well as w 2 . It could be originating from a larger initial eccentricity, 
although there should be a bias from selecting the larger u 2 events because the radii 
modulations also depend on the anisotropy in the momentum space ISO- 
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Figure 2: (Left and center) Azimuthal angle dependence of R% and R% relative 
to \I/ 2 . (Right) 2 nd -order azimuthal oscillation on R? s and R? a . In both panels, 
results with and without higher Q 2 selection are shown. 
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3 Study on the event twist effect 

The data of Pb+Pb 2.76 TeV collisions simulated using a Multi Phase Transport 
Model (AMPT, v2.25 with string melting) were used for the study on the twist effect 
of the particle-emitting source in the final state. The impact parameter was fixed 
to 8 fm. For the HBT study, the interference effect between two identical particles, 
1 + cos(Ar • Ap), was calculated and weighted to the relative pair momentum 
distributions. Then the correlation functions were reconstructed by taking a ratio 
of the distributions with and without the weight. Also, all charged pions were 
allowed to make a pair with each other including 7r + 7r~ to increase the statistics 
(the consistency between results for positive and negative pairs was checked). The 
event plane was determined using particles in 4 < |ry| < 6, where particles were 
divided into two sub-groups. A set of forward and backward event planes ('J'f,'!'®) 
were used for the event cut which requires finite difference between \E'| and tpf as a 
event twist selection, whereas the other set of event planes was used for a reference 
angle of azimuthal HBT measurements. The effect of the event plane resolution 
was not taken into account in this study, assuming the resolution hardly affects the 
phase of the radii oscillations but the magnitude of the oscillations. 

Figure [3] shows R^, R°, and R° s as a function of azimuthal pair angle relative 
to the (A 0 ) for four 77 regions, where ('F® — ^2) > 0-6 was required. The 
oscillations of the HBT radii measured in 77 > 0 are shifted to the direction of 
A tp < 0, which is the direction of 4/!, in the current event cut. This phase shift can 
be understood to be a possible twist effect in the final source distribution. The radii 
oscillations measured in 77 < 0 have no phase shift because the negative 77 is closer to 
the backward angle used for a reference of the HBT measurement and therefore the 
twist effect is supposed to be small. Figure[4]shows all the extracted radii parameters 
measured with respect to 4/^ and 4/f 1 with the event cut of (4/^ — 4/f) > 0.6. The 
phase shift can be seen in the transverse radii (not in the longitudinal radii including 
the cross terms) when comparing the results for 4/^ and ’I'.f in the same 77 region. 
It indicates the twist of the event plane angles in the backward and forward angles. 
These oscillations were fitted with the following functions: 

^(A0) = R l,o + 2 R l,2 COS (2A0 + a) (for 77 = o,s), (2) 




Figure 3 : Azimuthal angle dependence of Rj?, R%, and R^ s relative to the 
backward 4> 2 with the event cut of ( 4 >f — d'f) > 0 . 6 , where the dashed lines 
show i?Q S = 0 and the dotted lines show A <f> = tt/2. 
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Figure 4 : Azimuthal angle dependence of the radii parameters relative to the 
backward and forward v f , 2 with the event cut of (tF® — \Ff) > 0.6, where the 
dashed lines show = 0. 


^(Ac/>) = R l,o + 2R l,z sin(2A^ + a) (for fj, = os), (3) 

to extract the magnitude of the phase shift (a). The phase shift parameter a 
obtained from the results with respect to and rFf is plotted as a function of 

V in Fig. 0 The a increases going from backward to forward angle in all cases. 
The variation of a in the ?? dependence is comparable to the difference between 
results relative to \Ff and at the same rj. The finite slope in the r] dependence 
of a indicates the twisted source in the final state due to the initial longitudinal 
fluctuations, as well as the twisted event plane and flow as discussed in Ref. dug. 
In this study, the effect of the event plane resolution is assumed to be negligible, 
but there should be an uncertainty derived from the event cut of ('F® — ’Ff) and 
their finite event plane resolutions. 

4 Summary 

We presented the results of HBT measurements using the event shape engineering 
for Au+Au collisions at yts NN = 200 GeV at the PHENIX experiment. We found 
that the higher Q 2 selection enhances the measured final source eccentricity as well 
as V 2 - Although the model comparison is needed to disentangle both spatial and 
dynamical effects on the HBT radii, this study would clarify the relation between 
the initial and final source eccentricities and constrain better the system dynamics. 
We also studied the event twist effect using AMPT model. When selecting 
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Figure 5 : Phase shift parameters a obtained from the RR°, and R os as a 
function of r). Results measured with respect to forward and backward event 
planes > ^2) are shown, with the event cut of — 'Pj) > 0.6. 


events with finite difference between forward and backward event plane angles, the 
oscillations of HBT radii are shifted in the phase and the phase shift increases with 
77 . The results indicate a possible twisted source in the final state preserving the 
initial twist due to the longitudinal fluctuations. This effect might be measured in 
experiments at RHIC and the LHC. 

Both techniques could be useful to probe and control initial fluctuations in 
transverse plane and longitudinal directions, as well as to study the response of the 
system to the space-time evolution. 
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